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Chapter 1 

Dissipative Kerr solitons in 
optical microresonators* 

Tobias Herr^’^ 

Michael L. Gorodetsky^’^ 

Tobias J. Kippenberg^ 


Abstract: This chapter describes the discovery and stable generation of tem¬ 
poral dissipative Kerr solitons in continuous-wave (CW) laser driven optical mi¬ 
croresonators. The experimental signatures as well as the temporal and spectral 
characteristics of this class of bright solitons are discussed. Moreover, analytical 
and numerical descriptions are presented that do not only reproduce qualitative 
features but can also be used to accurately model and predict the characteris¬ 
tics of experimental systems. Particular emphasis lies on temporal dissipative 
Kerr solitons with regard to optical frequency comb generation where they are 
of particular importance. Here, one example is spectral broadening and self- 
referencing enabled by the ultra-short pulsed nature of the solitons. Another 
example is dissipative Kerr soliton formation in integrated on-chip microres¬ 
onators where the emission of a dispersive wave allows for the direct generation 
of unprecedentedly broadband and coherent soliton spectra with smooth spec¬ 
tral envelope. 


1.1 Introduction to optical microresonator 
Kerr-frequency combs 

Ultra high quality-factor (Q) optical whispering-gallery mode resonators have 
long been considered promising for their nonlinear optical properties since their 
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discovery by Braginsky and co-workers in 1989 [T]. Early works included the ob¬ 
servation of bistability, thermal nonlinearity [5], Kerr-nonlinearity at cryogenic 
temperature [3], stimulated Raman- [3] and Brillouin-scattering[S], electro-optic 
effect [6l [7j, as well as second[8] and third[9] harmonic generation. In partic¬ 
ular, the third order Kerr-nonlinearity can give rise to parametric oscillations 
in an optical microresonator. In this four-photon process two pump photons 
are converted to one signal and one idler photon (i. e. degenerate four-wave 
mixing) where the photonic energy is conserved. These parametric oscillations 
(sometimes called hyper-parametric to distinguish them from three-photon os¬ 
cillations associated with the second order optical nonlinearity), despite being 
well-known in nonlinear optics for decades |10|. were observed in ultra high-Q 
toroidal fused silica and crystalline resonators only in 2004 [ni [ID- These studies 
demonstrated the dramatic reduction of threshold power for nonlinear optical 
oscillation owing to the 1/Q^ scaling of the threshold[T]. 



b Frequency Domain 

|E(v)| 



Figure l.I: Time and frequency domain picture of mode-locked laser based fre¬ 
quency combs. A periodic train of pulses with a pulse repetition rate /rep (panel 
a) corresponds to a comb spectrum of equidistant lines in the frequency domain 
(panel b). The line spacing is given by /rep. The offset /g of the frequency comb 
spectrum relates to the carrier-envelope phase shift A(f>Q between two consecu¬ 
tive pulses via /g = /rep • A^o/(27r). The two parameters /rep and /g fully define 
all comb frequencies /n = n • /rep + Jo¬ 
in 2007 it was discovered [T3] that Kerr-nonlinear optical microresonators can 
give rise to optical frequency comb[T31 [H] generation. The surprising cascade 
and mode proliferation associated with four-wave mixing (FWM) was demon¬ 
strated to lead to a broadband series of optical lines equidistantly spaced in the 
optical frequency domain, that is an optical frequency comb. The discovery of 
these so-called microresonator Kerr-combs established optical microresonators 
as tools for precision frequency metrology, a connection not made before. Such 
frequency comb spectra are conventionally derived from mode-locked femto¬ 
second lasers. In this way, Kerr-combs broke with the conventional dogma in 
the frequency metrology community, that an optical frequency comb requires 
a mode-locked pulsed laser source. Different from conventional mode-locked 
laser based combs, regular microresonator Kerr-combs are not pulsed in the 
time domain. By comparing the comb spectrum from a microresonator with 
a conventional mode-locked fiber laser comb, the equidistance of the comb 
lines could be proven at the level of 1 part in lO^’^. Since this work the 
field of microresonator Kerr-combs has increased substantially. In the subse¬ 
quent years, Kerr-combs have been demonstrated in a variety of platforms, 
including crystalline resonators[T5], CMOS compatible platforms such as sili- 
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con nitride (Si3N4)[T71 [THl [TO] , Hydex glass [ID], as well as aluminum nitride 
|21| or diamond |22|. In addition to allowing for miniaturization and chip-scale 
integration of frequency comb oscillators, microresonator based Kerr-combs en¬ 
able to attain wide comb line spacings in the technologically relevant 10 — 100 
GHz range[23l |24|, which is not easily accessible using mode-locked laser fre¬ 
quency combs. Such large mode spacing is of particular interest in coherent 
telecommunication|25j . spectrometer calibration for astronomy uniiizj, low noise 
microwave |28l IDD] and arbitrary optical waveform generation |30j. Optical tem¬ 
poral dissipative cavity Kerr solitons, predicted theoretically jS!] [32] and first 
realized experimentally in optical hbers[33]. have played a decisive role in al¬ 
lowing to overcome noise processes associated with comb formation in microres¬ 
onators. As detailed later, temporal dissipative Kerr solitons (DKS) provide a 
deterministic way to generate broadband, coherent, and spectrally smooth opti¬ 
cal frequency combs, which are in addition amenable to an analytical theory and 
numerical simulation. These developments, along with the experimental genera¬ 
tion of single soliton states in crystalline resonators [341135] and integrated Si 3 N 4 
resonators |36]. make microresonators suitable for precision frequency metrology 
and many related applications. Owing to the pulsed nature of the solitons (in 
contrast to the earlier Kerr-combs), microresonators can now provide equiva¬ 
lent counterparts to mode-locked laser regarding the time domain properties. 
This allows using methods of spectral transfer and broadening as well as self- 
referencing techniques [SZ1IM1I3S] that have been developed for mode locked 
laser-systems. In this way DKS in microresonators provide a route to syn¬ 
thesize absolute optical frequencies from an RF or microwave signal and the 
opportunity to use microresonators for counting the cycles of light. 


1.2 Resonator platforms 

Parametric frequency conversion in Kerr-nonlinear microresonators, first ob¬ 
served in toroidal fused silica and crystalline microresonators, has now been ob¬ 
served in a wide variety of geometries and platforms, including microspheres [JD] 
m] and integrated, planar ring-type resonators, such as those based on Si 3 N 4 [ 32 ]- 
The work presented in this chapter on Kerr soliton generation has been carried 
out with crystalline magnesium fluoride MgF 2 resonators and Si 3 N 4 integrated 
microrings, that are briefly described in the next sections. 

1.2.1 Ultra high-Q MgF 2 Crystalline microresonators 

The first microresonator platform in which temporal DKS have been observed 
experimentally]^ were crystalline resonators made from magnesium fluoride(MgF 2 ). 
Crystalline optical resonators were first introduced in [43] and further develop¬ 
ment demonstrated that polishing of crystalline materials can lead to exception¬ 
ally high Q-factors (exceeding 10^^ [HJIIS]). This method therefore extended 
the ultra high-Q attained with surface reflow methods [TJHD] to crystalline ma¬ 
terials. Figure |1.2| shows an ultra high-Q crystalline resonator made from a 
MgF 2 cylindrical preform, which is prestructured with protrusions using pre¬ 
cision diamond turning. After prestructuring a series of polishing steps using 
diamond films and slurries are used to attain a low roughness surface finish and 
ultra high-Q optical modes. To excite the whispering-gallery modes evanescent 
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Figure 1.2: Microresonator platforms, (a) Diamond turned magnesium fluoride 
resonator containing two protrusions that confine high-Q optical whispering- 
gallery modes, (b) Manually shaped and polished resonator, (c) Polishing of a 
diamond turned preform, (d) Scanning electron micrograph of a silicon nitride 
microresonator before applying the fused silica cladding. 


optical coupling via a tapered optical fiber can be used [471 US]- This is possible 
as the refractive index of MgF 2 is lower than that of fused silica fibers. The use 
of thin tapered fiber for coupling is feasible despite the large volume of the em¬ 
ployed (typically mm scale) resonators, due to the exceptionally high-Q, which 
reduces the necessary coupling rates required for critical coupling. Typically the 
resonance width is in the range from 50 — 500 kHz corresponding to Q values of 
10® to several 10®. The nonlinear frequency conversion threshold is reached at 
optical powers below 1 mW. For soliton generation pump powers of the order 
of 10 mW are typical. 

1.2.2 Integrated photonic chip microring resonators 

From a future application perspective microresonator platforms amenable to 
wafer-scale processing can be advantageous. A key challenge in this context is 
attaining sufficiently high-Q for efficient nonlinear parametric frequency conver¬ 
sion. One platform that is particularly well suited is based on silicon nitride 
(Si 3 N 4 ) waveguides embedded in fused silica, is a material already used in the 
micro-electronic industry (part of the CMOS process) and is well suited for 
integrated photonic waveguides. Its refractive index of n « 2 enables tight op¬ 
tical confinement waveguides and microring resonators and the high bandgap 
(~ 3 eV) mitigates two (or higher order) photon absorption in the telecommu¬ 
nication band. While the achieved Q-factors (~10®) in integrated devices are 
still many orders of magnitude below that of polished crystalline resonators , 
higher Kerr-nonlinearity and tighter optical confinement, lead to high effective 
nonlinearities that are almost three orders of magnitude larger than in crys¬ 
talline resonators. This has enabled parametric oscillations to be accessible 
with pump powers of ^lOOmW as demonstrated first in Refs. [13149]. Sub¬ 
sequent work demonstrated low phase noise comb ODeration[50j via nonlinear 
synchronization, as detailed later. In addition Si 3 N 4 microresonators were the 
first microresonator platform that allowed for the observation of soliton induced 
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Cherenkov radiation [3T] (or dispersive wave emission), i.e. the dynamics of soli- 
tons in the presence of higher order dispersion]^. Fabrication of the resonators 
proceeds by lithography, etching and a final encapsulation technique in which 
the waveguides are clad with fused silica. The image of a final (but not yet clad 
resonator) is shown in Figure 1.2 Coupling to and from the chip is achieved 


with inverse tapered waveguides, adiabatically converting the mode spot diam¬ 
eter from a tight nanophotonic waveguide mode to a large area mode that can 
be excited with a lensed single-mode optical fiber. 


1.3 Physics of the Kerr-comb formation process 


In a simplified picture, the Kerr-comb formation in microresonators starts with 
an initial degenerate FWM process, generating a pair of sidebands symmetri¬ 
cally spaced in frequency around the pump laser in the resonances adjacent to 
the pumped resonance. In this process the free spectral range (FSR) defines the 
frequency spacing of the sidebands from the pump laser. A subsequent cascade 
of non-degenerate four-wave mixing processes generates a large series of side¬ 
bands in frequency steps corresponding to the initially defined line spacing. The 
bandwidth of the comb is limited by the dispersion of the resonator that causes 
a mismatch between the optical sidebands and the resonance frequencies. This 
mismatch increases with the spectral distance to the pump wavelength. While 
it explains well the early results of Kerr-frequency comb generation with typi¬ 
cally THz-line spacing and moderate bandwidth, this picture had to be refined 
when attempting broader optical bandwidth and narrower line spacing (below 50 
GHz) in toroidal [53], crystalline [16] and integrated Si 3 N 4 resonators [TTllMl [55] . 
Here unexpected noise phenomena were observed in the form of broad optical 
lines, high amplitude noise and loss of coherence. These phenomena appeared 
to be independent of the chosen microresonator platform. Indeed, the origin 
of this noise can be explained by the universal properties of the comb forma¬ 
tion mechanism [50]. While in early microresonator comb experiments, due to 
the wide line spacing, the mode proliferation occurred indeed on adjacent res¬ 
onators modes, the latter is typically not the case in resonators with narrow 
FSR. Here the first pair of primary sidebands generated by the initial FWM 
process can occur in resonances that are widely separated from the pump (by 
multiple FSR). The process is depicted schematically in Figure 1.3 Using the 
nonlinear coupled mode equations, as introduced in the next section, one can 
show that the mode spacing between the primary modes is essentially defined 
by the (normalized) ratio of second order dispersion parameter and resonance 
width. If the resonator has a comparatively narrow FSR, then the initial side¬ 
bands are separated by a frequency interval A that encompasses a high number 
of FSR intervals. Further four-wave mixing can fill up the unpopulated reso¬ 
nances by formation of secondary sidebands separated by a smaller frequency 
spacing S. In other words, sub combs with d-line spacing form around the A- 
spaced primary lines. As, however, the primary spacing A is not guaranteed 
to be an integer multiple of the secondary spacing 5, the generated lines can 
not generally form a consistent comb, which manifests itself in the previously 
mentioned noise. This hypothesis of comb formation has been experimentally 
tested in |50j using a fiber laser frequency comb as a reference to reconstruct 
the frequency components of the Kerr-comb. Another key finding of this work 
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is that in contrast to the early model of comb formation, more than one single 
comb line can occupy a given resonance (cf. Figure 1.3 and Figure [l^. 
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Figure 1.3: Universal Kerr-comb formation processes, (a) Formation of pri¬ 
mary sidebands (b) Formation of subcombs (c) Overlap between inconsistent 
subcombs can lead to multiple lines per cavity resonance and explains noise 
phenomena in Kerr-combs. 

Despite these noise processes, regimes have been found where broadband 
coherent Kerr-combs can be generated. The first observed transition to such a 
regime of low phase noise was presented in m- Here the subcombs defined by A 
and S can be synchronized by changing the pump laser frequency and power so 
that A and S become commensurate. Once close enough to the commensurate 
state, nonlinear synchronization and locking sets in and creates a phase coherent 
optical frequency comb [50] (albeit not pulsed in the time domain). Such syn¬ 
chronization has been experimentally observed and studied by reconstruction of 
the frequencies of the Kerr-comb lines as shown in Figur ^TS] Once the Kerr- 
comb falls into such a synchronized state it remains stable against fluctuations 
in the pump laser parameters and can be used for further experiments. 

Later experiments have revealed that the nonlinear locking process can 
be captured by the Adler equation, that is known from injection locking of 
lasersjSSj. Kerr combs driven into a low-phase noise, coherent state have been 
used in applications, such as coherent communications [25] . where they are used 
to provide a multitude of data channels from a single laser, enabling terabit 
per second data-communication. Yet, despite the promising prospects a cen¬ 
tral missing element has been ways to reliably achieve such low phase noise 
combs, with smooth spectral envelope. A surprising discovery has been made 
in 2012(531311: when it was demonstrated that the generated subcombs can not 
only undergo nonlinear synchronization, but moreover can seed the formation of 
temporal DKS in optical microresonators. This discovery provides not only an 
opportunity to study DKS in microresonators, but moreover a powerful tool to 
generate low phase noise combs with broad spectral bandwidth that can be used 
in a variety of applications. Before turning to soliton formation, we review first 
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Figure 1.4: Kerr-comb formation and noise, (a) Formation of the optical spec¬ 
trum as the laser is tuned into resonance and the intracavity power increases, 
(b) Comb-line spacing measured as the radio-frequency (RF) beatnote between 
neighboring comb lines. Multiple and broad beatnotes indicate multiple and 
inconsistent line spacings present in the comb spectrum. 


the basic formalism describing Kerr-nonlinearity driven parametric oscillations. 


1.3.1 Nonlinear coupled mode equations 

The theoretical analysis of parametric oscillations, frequency comb and soliton 
formation is a complex nonlinear mathematical problem for which several ap¬ 
proaches have been proposed [niiMiiiiiiMiiiniEi]- Besides description in the 
time domain the modal expansion approach in the frequency domain [621l63j has 
proven particularly useful in the context of microresonators as it allows easily 
to take into account particularities in the mode structure. The system of equa¬ 
tions describing the dynamics of each optical mode Am in the slowly varying 
envelope approximation can be obtained from the nonlinear wave equation in 
conjunction with the quantum Langevin equations: 


dt 


_Km . . 

2 -^m ~ ^m,mo 




huJn 


ig ■ 




jyr. 


Am' Am" Am'" C 


m' ,m" 


4” (^) 4“ \/(^) ■ 


( 1 . 1 ) 


Here Km denotes the cavity decay rate of mode m with eigenfrequency ujm, ^ 
is the reduced Planck constant. The cavity decay rate is composed of intrinsic 
(Ki^m) and external coupling (Kex.m) loss rates such that the total loss rate is 
Km = Ki,m 4- Kex,m = , where Qm Is the loaded quality-factor of mode m. 

One mode (mo, wo) is driven by a classical field with optical power Pin and 
frequency ujp. The terms Ssi^m and Ssex,m designate quantum Langevin input 
noises with commutation relations [5sj(t), i5sj(t')] = S{t — t') and with zero 
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Figure 1.5: Subcomb synchronization in a Si 3 N 4 optical microresonator achieved 
by changing the pump laser wavelength.. 


mean {{Ssj(t)) = 0). The summation over the nonlinear interaction terms is 
made over all azimuthal indices for the nearly equidistant set of fundamental 
whispering-gallery eigenmodes, satisfying the following conservation relation: 
m' + m” — ml" = m, where is the degeneracy factor: I?™ = 2 in case 

of cross-phase modulation and non-degenerate four-wave mixing, and I?"* ™ = 1 
for self-phase modulation and hyperparametric generation. We neglect processes 
where one of the frequencies is significantly smaller (ujm = ^m' — ^m" — ) 

or larger than other frequencies (wm = + Wm'") as it is the case for 

third harmonic generation (when in the degenerate case 2?^™ = 1/3). The 
intermodal coupling factor is given by 

™ “ /|eo|W x(3)(a;o = a;o+u;o-cuo) ’ 

where is vector electric field distribution of the mode m and describes 
the third order Kerr-nonlinearity. This factor is analytically different from that 
obtained in [611 Ell |63] The effective nonlinear volume of the modes is given 
by: 


^ (/ \em?dVf 

/|e™|W 


(1.3) 


The latter may be found from the cavity eigenmodes’ field distribution. The 
modes are normalized such that corresponds to the number of photons 

in a given cavity mode. The nonlinear coupling constant g (that is the Kerr 
frequency shift per photon) is defined as in [58] : 


_ (‘^o) _ H.ujnCn2(ujo) /, 

y 4£on"‘(i.iJo)Vb \ ) 

^The authors of I61ll62|[63l overlooked in the derivation that the nonlinear polarization is 
time dependent and can not be simply taken out from the time derivative in wave equation 
1641 . This omission leads to unphysical asymmetry in A over mode indices and loss of pre¬ 
cise degeneracies important for analytical considerations, but usually plays negligible role in 
numerical simulations of narrow combs when A ~ 1. 
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Here eg denotes the dielectric constant and n and n 2 are the refractive non¬ 
linear optical indices. A seeding via input noises is important for numerical 
simulations to initiate the comb, however, it is not required in further analyti¬ 
cal considerations as their means vanish. 

For the calculation of the field distribution and eigenfrequencies in crys¬ 
talline resonators we do not use the modes of a sphere as proposed in |62[ 163) 
but asymptotic solutions obtained for a spheroid [551|M|[S7]. producing better 
approximations for the wide range of geometries, or eigenfrequencies obtained 
via CDMSOL finite element modeling. From these approximations as well as from 
Sellmeier’s equation for the resonator material one can determine Vm and eigen¬ 
frequencies ujrn to find second and higher order dispersions. The latter allows 
can be approximated by 

where ^ = m — toq. The parameter D 2 is related to the second order group 
velocity dispersion (GVD) (32 {Pj = via D 2 = —^DfP 2 , and D 3 relates to 
the third order dispersion via D 3 = —^DfP 3 + 3^Z?i/3| « — 

It is convenient to remove the explicit time dependence in the nonlinear 
equations with the substitution where I?i/(27r)is the 

mode spacing of the comb at the pump wavelength. Moreover, the pump term is 
written as F = ^ ■ In this new basis the oscillations of each optical mode 

are considered not around eigenfrequencies but around nearest frequencies on 
an equidistant Di spaced-grid: 

~ DifJ, — cOp) + -^]Ap + 5fj,oF 

tg (1.6) 



If L pairs of sidebands and the pump are considered, then the total number of 
nonlinear terms in all 2L -|-1 equations is \{L + 1)(8L^ -I- 7L + 3). 

It is possible to get some useful conservation laws for this set of equations. 
Calculating the evolution of the total number of photons in the comb ^ = 




dN^ 
L dt 


with 


dN^ 

dt 


each other: 


-K" 


dAp 

dt 


— A 

— jj. -Ai- 


dA* ' " 

^ , -^Ap_ and assuming symmetry: ^ ^ -|- 

we find that most of the terms in the sum cancel 


dN 

dt 


L 

+ ^ KpNp = 

^— — L 


(FAS + F*Ao), 


(1.7) 


which in case of = k simplifies to 

dN 

-H kTV = 2|F||Ao| cos(()o, (1-8) 

where pQ is the phase difference between F and Aq. This master equation for 
the whole optical spectrum may be useful for example for the analysis of the 
steady states. 
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In a similar way one may find: 

^ ^ ^ ^ /r (k+^xA^+m - ^ (1-9) 

ji— — L ^ ^ fi——L 

which in the case of steady state means that the ’center of mass’ (i.e. the fre¬ 
quency center of photonic energy) of the spectrum is conserved and coincides 
with the pump frequency as each four-wave mixing process follows this conser¬ 
vation law. 

In the following analysis we use further simplifications, assuming equal pa¬ 
rameters for all comb lines A™ ™ ™ =1) and renumbering comb 

lines with index /i, starting from the central driven mode rriQ, i.e. /i = m — toq. 
In this case it is possible to write the equations in a dimensionless way, as 
suggested in [HH] : 

fl/x = 

f = ^/8g/K^\Fp\, 

C;x = 2{ujp,- flDi-UJp)/KKiCo+D 2 fl‘^/K, 

T = Kt/2, (1.10) 

doj 

~[1 + ^ ^ F Sop^f. (1-11) 

/x',/x" 

In this form all frequencies, detunings and magnitudes are measured in units 
of cavity resonance width so that = 1 corresponds to the nonlinear mode 
pulling of one cavity resonance width (thresholds for both single mode bistability 
and degenerate oscillations). 

1.3.2 Degenerate hyperparametric oscillations 

An analytic solution of the complete nonlinear system with many comb lines 
in general case is hardly possible. It is possible, however, to resolve equations 
for an arbitrary set of modes if only degenerate processes caused by the pump 
as well as nonlinear mode pulling effects corresponding to cross- and self-phase 
modulation are taken into account. These effects are dominating as long as the 
energy in each of the optical lines is small compared to the intracavity pump 
power and when due to strong second order dispersion the resonator modes are 
far enough from equidistance lines to suppress non-degenerate four-wave mixing. 

The following set of reduced equations qualitatively describes the processes 
of initial comb formation, including widely spaced primary lines, switching be¬ 
tween initially excited modes, chaotic behavior and noisy RF signals with many 
frequencies: 

Bn 

= ~[1 + *Co]ao + i{‘2N — |aoP)ao -I- 2i o.pO,-pO,Q + /> 

So, ~ 

-^ = -[^ + iC,p\ap + i{‘2-N-\ap\'^)ap + ia*_p^al, ( 1 . 12 ) 

where according to the normalization N = 2gN/n. The reduced model was com¬ 
pared with the full numeric model integrated using the Runge-Kutta scheme. 
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It shows good agreement concerning the initially generated sidebands and qual¬ 
itative agreement of the comb in the noisy, chaotic regime. 

With the substitution = 2(/)o — cj)^ — cj)-^ it is possible to 

search for steady state solutions with = 0 and = const. One may 

obtain a a steady state from the resulting equations when (/jq = 0 and = 0, so 
that = —(p-fi and Soj^^ = —6u}-^ = /(in) -|-.... Splitting the equations 

(1.12) into real and imaginary parts we get: 


ao(l -l-2^a^Q!_^sin'0;i) = f cos4>o, 

L 

(Co - 21V - 2 ^ O/iO-zx cos + al)aQ = / sin < 

sintjjf, = 0, 

Cm«/.* - (2^ - oil)af, - cos V'/x = 0 , 


(1.13) 


where 


C/x = (C/x - 20 ;x/k)- 


(1.14) 


From the third equation we find that and tpij, = ip = arcsin(l/aQ). 

Substituting this into the first and second equation of (1.13) we obtain: 


N = fao cos (po, 

[Co -2N + al + {N - ao)\/l - 1/Q!o]«o = fsincpo- (1.15) 


The first equation of (1.15) also follows from the master equation (1.8). Adding 


both equations squared we get the system: 


[Co ~ N(2 + cos Ip) + cxq{1 cos ip)]"^(X q /Og — 
^ , „2 ' ~ 


a“ - 2iV = CrxZ = CtQ COS Ip. 


(1.16) 


It follows that CnZi which determines the detuning of the new nonlinear 
effective eigenfrequency from the frequency of oscillations, does not depend on 
fjL. This means that in the steady state an equidistant comb may be formed if the 
dispersion of each mode is compensated by nonlinear pulling with appropriate 
amplitude distribution of each excited comb line. This nonlinear dispersion 
compensation is possible due to the fact that self phase modulation is two times 
smaller than cross phase modulation and is an essential mechanism for comb 
formation uni EH]. 


The last equations of (1.16) may be summed over all excited lines to get: 


N = 


Cs - ^0 + 2Av^ao - 1 


4L- 1 


(1.17) 


where Cs = C = (2A + l)Co + d 2 (where d 2 = D 2 /k). If excited 

comb lines are filling all resonator modes in some interval without omissions 
the summation may be done analytically. If substituted in the first equation of 
(1.16) we finally can get an equation for the pump mode magnitude a§ ready 
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for numerical analysis. Substituting /3 = Oq — a/oq — 1, then trans¬ 

forms the final equation to a rational 6-th order equation. Further analysis of 
these equations may reveal the properties of the low-noise low-populated combs 
observable in experiments. 

1.3.3 Primary sidebands 

A particular case of the reduced model is a three mode system for which all the 
nonlinear terms are in fact considered, i.e. the model is complete. This system 
can describe the threshold for primary sideband hyperparametric oscillations at 
a /i times FSR frequency separation from the pump that can be softly excited 
when the pump is gradually tuned into resonance starting from infinite detuning 
{ujo < ujp for I ?2 > 0). 

The three mode system and its stability was analyzed in detail using numer¬ 
ical analysis of the nonlinear equations |58j or stability analysis of the linearized 
equations [52]. Both approaches lead to the same equation following from the 
reduced model when L = 1: 

From the nonlinear resonance condition and condition for the first sideband 
to emerge obtain: 

{C,-oclfo.l + al = f\ 

Co + d 2 ^? = 2al - (1.18) 

It follows the equation determining the possibility to excite the mode /r: 


Co = ao - \IPIal - 1, 

\l- 1 - d2lith + «o - \/ao- 1 = 0- 
or using the smallest possible Og = 1 when the radical is real: 


/^th 


7(1 +V/2-i)M 


Alternatively, this can be rewritten in physical terms |50j as: 


Mth — 



(1.19) 


( 1 . 20 ) 


( 1 . 21 ) 


The threshold power Pth corresponds to |ao,th| = 1- The minimum obtainable 
for Pin = Pth is given by: 


Mth = 



( 1 . 22 ) 


and depends only on the ratio of cavity decay rate k to second order cavity 
dispersion P> 2 - 

Note that the threshold oo.th = !> determining the first sideband to appear 
when tuning the pump into resonance at a given power (also found in [62] 1. 
however, does not correspond to a minimum of input power when hyperpara¬ 
metric generation may start. This minimum, numerically calculated in [58] may 
be found explicitly as Oin.min = 2/^3. 
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1.4 Dissipative Kerr solitons in optical microres¬ 
onators 

The formation of stable temporal DKS circulating in a continuous-wave (CW) 
driven microresonator with Kerr-nonlinearity is a remarkable form of optical 
self-organization. It requires, similarly to the Kerr-combs described above, a 
medium with focusing third order nonlinearity and anomalous group velocity 
dispersion. Fundamentally different from the above described Kerr-combs, the 
formation of temporal DKS implies a stable pulsed time domain waveform where 
the pulse properties are fully determined by the experimental parameters. The 
ultra-short soliton pulse duration, typically in the femto-second regime, implies 
broadband optical frequency comb spectra that show the same wide line spacing 
that is known from the previous Kerr-combs. Importantly, however, the soliton 
based comb spectra do not suffer from the inconsistent subcomb formation and 
associated noise processes that are often observed in Kerr-combs. Prior to the 
observation in a microresonator, temporal DKS had already been observed in a 
several hundred meter long, CW driven fiber cavity, where additional short laser 
pulses had been used to stimulate the soliton formation |33j . In this context it 
was speculated that DKS could be formed in a microresonator as well (Note: 
the possibility of solitons circulating in a dielectric sphere was already proposed 
in Ref. [69]). Moreover, Refs. (TO] [71] pointed out that the theoretical formal¬ 
ism so far used by researchers to describe microresonators and fiber-cavities, 
i.e. the non-linear coupled mode equations and the Lugiato-Lefever equation 
m respectively, can be mapped onto each other. From a historical perspec¬ 
tive, it is thus interesting to note that after the first reports on microresonator 
Kerr-combs in 2007[T2], it took five years to enter the regime of Kerr solitons in 
a microresonator. One reasons for this is the particular thermal behavior of a 


microresonator system that had long masked this regime (cf. Section 1.6 1 . More¬ 


over, the approach of stimulating solitons using laser pulses that is successfully 
used with long fiber cavities can not easily be transferred to microresonators 
due to their much higher finesse (i.e. difficulty of coupling the pulse to into the 
resonator). Before presenting experimental observations, the following section 
will establish the theoretical formalism that allows describing dissipative Kerr 
soliton in a microresonator. 


1.4.1 Analytical theory of dissipative Kerr solitons 

To describe the internal field in a nonlinear microresonator the following mas¬ 
ter equation may be used when third and higher order dispersion terms are 
neglected [SUl 170] : 


dA .1 d^A 


zg\A\^A 


+ f(wo - UJp 


A 


/ 


nrjPjn 

hujQ 


(1.23) 


Here A{(j),t) = is the slowly varying field amplitude and 

(j) € [0, 27r] is the angular coordinate inside the resonator with periodic bound¬ 
ary conditions, rj = Kgxl^ is the coupling efficiency with critical coupling corre¬ 
sponding to rj = 1/2. Transforming eq. (1.23) to its dimensionless form gives: 
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Figure 1.6: Crystalline resonator for dissipative Kerr soliton generation, (a) 
Crystalline MgF 2 resonator and simulated mode profile, (b) Measured anoma¬ 
lous group velocity dispersion for the resonator shown in panel a. The dispersion 
of a resonator can be quantified in terms of the deviation of its resonance fre¬ 
quencies w^from an equidistant frequency grid loq + fJ-Di, where Di/{2tt) is the 
FSR at the pump wavelength. An anomalous group velocity dispersion corre¬ 
sponds to a parabolic curve as shown in panel b. 


Here 4'(r, (/>) = X) is the waveform, 6 = is the dimensionless 

longitudinal coordinate, and d 2 = D 2 /H is the dimensionless dispersion. Equa¬ 
tion (1.24) is identical to the Lugiato-Lefever equation where a transversal 
coordinate is used instead of a longitudinal one in our case. It is important to 
note here, that though we used throughout the chapter the term temporal to 
differentiate DCS from better known spacial cavity solitons in nonlinear Fabry- 
Perot etalons, this term may also be confusing. For temporal dissipative solitons 
in fibers and fiber loop resonators in equations time is usually considered as a 
“fast” variable, describing the form of the pulse, while the coordinate along prop¬ 
agation is describing slow variations. In the case of WGM microresonators the 
opposite mapping is more convenient with ’’slow” time variation averaged over 
many roundtrips and fast mapped azimuthal angular coordinate. The equation 
1.24| may also be considered as a driven, damped Nonlinear Schrodinger Equa¬ 
tion (NTS) [31]. Analytical solution of the NTS in the form of bright soliton 
is known only in case of zero dissipation m- An approximate solution for 
the dissipative soliton may be found using Lagrangian perturbation approach 
[Zllig. Eirst we introduce Lagrangian density so that the variation of it over 
dt* (Landau-Euler equation) leads to an unperturbed NLS equation: 
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Now taking into account the perturbed equation (1.24) we introduce the dissi¬ 
pative function: 
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The equations that need to be solved are the following: 


dL d dL 
dvi dr dfi 
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(1.28) 


where ri are different possibly time dependent collective coordinates. Using the 
ansatz of a stationary = 0) soliton i?e®‘^“sech(i30) (exact for unperturbed 
case when B = with n = B and r 2 = (po [32] we get: 

L = -2B^ + ^B^ -2BCo- ( 1 . 29 ) 

This yields: 

dB 

—— = — 25 + tt/cos 7 ( 1 . 30 ) 

dr 

(1.31) 

This leads to the stationary parameters of the soliton attractor [32] : 


B = x/2Co, 

(1.32) 

25 PWo 

cospo = = 

TT/ TT/ 

(1.33) 
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(1.34) 


An approximate single soliton solution which accounts for a flat CW background 
is then given by 


~ + Se*‘^“sech(S6*), 


(1.35) 


The constant CW background dtp may be found by inserting vko into eq. (1.24) 
as the lowest branch [ST] of the solution of 


(|«'op-Co+*)^o = */, (1.36) 

which, eventually, results for Co > a/S (bistability criterion) and large enough 
detunings p < ^Co(Co + 9) in: 
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(1.37) 


Extending (1.35) to the case of multiple solitons inside the resonator gives 

'ko + +*^2Co-^) Ef=i sech(^(C>- <(.,)). (1.38) 


This ansatz automatically translates to discrete possibles states (steps) in non¬ 
linear resonant response. These steps are an important signature of soliton Kerr 
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combs and are observed in the experiments (see below). The height of steps in 
the intracavity power can be found by averaging the waveform amplitude (eq. 
1.38) squared over one roundtrip for different numbers K of solitons: 


1^1 


= 




It: Jo 
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1 + 




(1.39) 


The Fourier transform {F) of a hyperbolic secant soliton is again a hyperbolic 
secant: 





-\/d2/2sech 




(1.40) 


Using the relation for the optical frequency uj = ujp+fjDi and the time t = 4>/Di 
the spectral envelopes and the soliton waveform can be rewritten: 


4'(w — ujp) = \/d 2 / 2 sech((w — ujp)/Auj) with Aw 


2Di 

TT V ^2’ 


(1.41) 


and 


'l'(t) = sech{t/At), with At 



(1.42) 


The minimal possible soliton duration can be found by using 
the above equation for At: 


(eq.p3l 
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(1.43) 


This equation can be recast in form of the group velocity dispersion P 2 = 
^^ 02 /Df, the nonlinear parameter 7 = 7 ;^^ (for simplicity we assume criti¬ 
cal coupling 7] = ll2 and on resonant pumping): 


At 


min 


A 


(1.44) 


where denotes the finesse of the cavity. Note that the values Aw and At 

need to be multiplied by a factor of 2 arccosh(-\/2) = 1.763 to yield the FWHM 
of the sech^-shaped power spectrum and pulse intensity, respectively. 

For the time bandwidth product (TBP) we find At ■ Aw = 2 / 7 r or, when 
considering the FWHM of spectral and temporal power (in units of Hz and s), 
TBP = 0.315. 

In the case of K multiple solitons inside the cavity there is a more structured 
optical spectrum 4'^, resulting from interference of single soliton spectra 4'j(/r), 
where the relative phases of these spectra are determined by the positions (pj of 
individual solitons: 4>if(/r) = 'I>(/r) . The line-to-line variations can 

be high, however, the overall averaged spectrum still follows the single soliton 
shape and is proportional to V^4>(/r). 
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1.5 Signatures of dissipative Kerr soliton forma¬ 
tion in crystalline resonators 

The generation of DKS in microresonators is associated with several character¬ 
istic signatures that can already be observed during a transient generation of 
the solitons. Here transient generation refers to the situation where the pump 
laser is scanned (with increasing wavelength) through the resonance of the mi- 
croresonator. The described signatures are particularly useful from a practical 
perspective as they allow for rapid sample characterization. Depending on the 
detuning during the laser scan the system can enter different nonlinear regimes 
including, parametric sideband generation with Turing patterns [76], Kerr-comb 
formation, breather solitons m and stable dissipative Kerr soliton generation. 
One signature, though not unambiguously connected to soliton formation, is a 
transition from high intensity noise and and high phase noise (related to Kerr- 
subcomb formation) to a low noise state. While scanning the laser, the noise can 
be monitored by e.g. recording the intensity noise in the generated spectrum or 
by detecting the beatnote signal between the generated optical lines as shown 
in Figure PTtIj. 
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Figure 1.7: Signatures of soliton formation, (a) A staircase-like step structure in 
the pump laser transmission indicates the formation of several temporal DKS in 
the microresonator. High intensity noise is apparent in the transmission signal 
prior to soliton formation, (b) The formation of the solitons is associated with 
a transition to low noise. This is evidenced here by the transition of a broad RF 
beatnote (cf. Figure [T t] ) to a single narrow-band signal. (CF: center frequency) 

A more unique and easily observable signature of soliton formation in opti¬ 
cal microresonators is the occurrence of discontinuous staircase-like steps in the 
transmitted pump laser power while the laser is continuously scanned across the 
optical resonance (cf. Figure ) [2lQ When operating a microresonator a 
lower pump laser transmission implies a higher average intracavity power. The 
emergence of the steps may be explained by referring to the cavity bistability 
curve that describes the power inside the cavity in dependence of the pump 
laser detuning. When approaching the cavity resonance from the blue detuned 
side (laser frequency is higher than the resonance frequency and correspond¬ 
ingly laser wavelength is shorter than the resonance wavelength) the intracav- 

^Discontinuous steps can also be observed (typically with a higher signal-to-noise ratio) in 
the power of the nonlinearly generated light. 
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ity power will increase. This increase of intracavity power will then, due to 
the Kerr-nonlinearity, effectively shift the resonance frequency towards longer 
wavelength. As the laser detuning is further decreased the intracavity power 
will steadily increase (upper branch on the bistability curve), until the highest 
possible power is reached at the point of effective zero-detuning. Beyond this 
point (now effectively red-detuned) the intracavity-power will steeply drop to a 
much lower value (lower branch on the bistability curve) as the Kerr-nonlinear 
resonance shift vanishes (cf. Figure 1.8 1 ). Exactly at this transition to red 
detuning temporal DKS can emerge from the noisy intracavity waveform (cf. 
Figure 1.9). In this case the formation of a high peak power soliton pulse will 
cause an additional Kerr-frequency shift such that the part of the pump light 
that sees the soliton pulse (i.e. copropagating with the soliton inside the cavity) 
follows now a second bistability curve. |1.8|3 . The result is a second power drop 
at a larger detuning when the soliton vanishes. 




Figure 1.8: Stability of soliton states: Considering only the Kerr-nonlinear res¬ 
onance shift the intracavity power can be described by bistability curves where 
the upper branch solution corresponds to high and the lower branch solution 
to low intracavity power. When tuning into the resonance with decreasing op¬ 
tical frequency (increasing wavelength) the intracavity power follows the upper 
branch of the Kerr-bistability curve. After the transition to a soliton state the 
major fraction of the pump light is described by the lower branch of the bista¬ 
bility curve. The fraction of the pump light that propagates with the soliton 
inside the microresonator experiences a larger phase shift and is effectively blue 
detuned on the upper branch of another bistability curve. The extend of the 
’soliton bistability curve’ towards longer wavelength depends on the peak power 
of the solitons (i.e. the maximal nonlinear phase shift), the relative height of 
the curve depends on the relative fraction of the pump light that is affected by 
the high peak power soliton. The overall intracavity power can be inferred by 
adding the bistability curves resulting in the black curve. 


The combined signal of power drops occurring at different laser detunings 
leads to the overall staircase-like step structure in the transmitted power. If 
a higher number of soliton pulses is generated, a series of multiple steps can 
be observed (during a rapid scan not all solitons disappear at the exact same 
detuning). The exact number of generated solitons differs from one laser scan 
to the other as the solitons emerge from a noisy state. It is interesting to note 
that the transition to effective red-detuning has not been considered in previous 


work as it is usually associated with only low intracavity power m- Section [T6 
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will describe how stable operation in the soliton regime can be achieved despite 
the steps in the intracavity power that can affect the resonator temperature. 
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Figure 1.9: Laser detuning and soliton formation. Panel (a) shows the transmit¬ 
ted power and a series of steps associated with multiple temporal Kerr solitons. 
The background shading indicates the laser detuning that is derived from a 
Pound-Drever-Hall (PDH) error signal shown in panel (b). The soliton forma¬ 
tion coincides with the transition to red detuning (indicated by a sign change 
in the PDH error signal). Note that the additional Kerr-frequency shift due to 
the high peak power soliton does not significantly impact the PDH error signal 
as its relative contribution is small. 


1.6 Laser tuning into the Kerr soliton states 


For stable operation of a microresonator system thermal self-locking [79l [80] 
is a widely used method. Here the resonator self-locks itself thermally to the 
pump laser frequency. Besides being a simple and convenient method, this al¬ 
lows for a fixed pump laser frequency during steady state operation (instead 
of a pump laser that follows drifts in the resonator). This is advantageous for 
frequency comb operation, where the pump laser constitutes one of the comb 
lines that should remain fixed in frequency. The mechanism of thermal self¬ 
locking is somewhat analogous to the Kerr-nonlinear bistability curve and gives 
rise to a so-called thermal triangle, that is known from the earliest reports 
of whispering-gallery mode microresonators [Tl|2]. In a resonator with positive 
thermo-refractive coefficient (as in MgF 2 or Si 3 N 4 ) thermal locking is achieved 
by blue detuned pumping. In this regime an increase of pump wavelength (or 
equivalently a decrease of resonance wavelength) implies a higher intracavity 
power and higher rate of absorptive heating. As a consequence of this heating 
and mostly due to the thermo-optic effect the resonance wavelength shifts to¬ 
wards longer wavelength thus stabilizing the relative detuning between pump 
laser and resonator. Exactly this behavior is present during a soliton state. 


This can be seen e.g. in Figure 1.7 where (within a particular soliton state) an 


increase in pump wavelength causes an increase in intracavity power (decrease 
of transmission). 

In order to stably reach the regime of thermal locking of a soliton state it is 
important that the transition to the soliton state does not induce a change of 
the resonator temperature (despite the fact that the intracavity power changes), 
which would directly destabilize the system. One solution to this problem is to 
tune the pump laser from a far-blue detuned position into the soliton state with 
a tuning speed chosen such, that the “integrated absorptive heating” over the 
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tuning time corresponds to the heating rate in the soliton state. In this case 
the temperature does not change significantly when the soliton state is reached. 
Thermal locking then ensures the stable operation of the soliton state (for many 
hours), without actuation on the pump laser frequency, as illustrated in Figure 
|1.10| It is important to note that the solitons emerge from intrinsic noise present 
in the cavity prior to reaching the soliton state. As a consequence the initial 
number of solitons generated via the laser tuning method varies randomly from 
shot to shot. However, the experimental conditions can often be adjusted such 
that a certain number of solitons (e.g. one soliton) is generated with high 
probability. The maximal number of soliton that can be generated via the laser 
tuning method is related to the dispersion of the cavity and grows with the 
applied pump power |34] . For typical parameters only a few, maximally up to 
approximately 10 or 20, solitons are generated. As the exact number depends 
on a random process the electronically controlled laser tuning method is usually 
repeated several times until it yields the desired number of solitons. Note that 
the laser tuning method does not require an 

additional external stimulation (e.g. a laser pulse), as it is the case in fiber 
cavities. The following section will characterize the solitons generated via the 
described laser tuning method. 



Figure 1.10: Generation of stable solitons via laser tuning, (a) Illustration of 
laser transmission, scan control voltage (corresponding to pump wavelength) 
and resonator temperature for slow, ideal and fast laser tuning speed. In the 
ideal case the resonator temperature does not change once the desired soliton 
state is reached and the laser scan can be stopped for stable operation, (b) 
Regular scan of the pump laser over a resonance showing a ’soliton step’, (c) 
The laser tuning method allows to stop the laser scan once the soliton is gen¬ 
erated. Once generated in this manner, the soliton circulates stably inside the 
microresonator. 


1.7 Simulating soliton formation in 
microresonators 


Dissipative temporal Kerr solitons in a microresonator can be simulated by in¬ 


tegrating the coupled mode equations (1.111 in time similar to work related to 
Kerr-combs by [niMiini]- We use an adaptive Runge-Kutta method of the 5-th 
order, which together with a transformed system of equations without trigono¬ 
metric functions and explicit time dependence allows achieving fast simulation. 
The simulation of hundreds of optical lines allows for accurate modeling of real¬ 
istic experimental systems. This includes in particular the beatnote between the 
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optical lines (i.e. the repetition rate), the envelope of the optical spectrum as 
well as instant amplitudes, phases and frequencies of all modes. The sum of the 
nonlinear mixing terms is efficiently calculated in the Fourier-domainjSJ. An 
additional self-steepening term is added when simulating broad combs in Si 3 N 4 
rings to better model the behavior of few-cycles pulses [82] ■ In simulations, sim¬ 
ilar to the experiment, solitons can be generated via tuning the pump laser into 
resonance. The number of generated solitons depends on the noisy intracavity 
waveform preceding the soliton state. To allow for deterministic simulations of 
a certain number of solitons, the respective number of sech-pulses can be seeded 
as an initial waveform while the simulation is started directly with a detuning 
that allows for stable solitons. It is worthwhile noting that the underlying cou¬ 
pled mode equations readily allow for the simulation of complex and irregular 
mode structures as often encountered in real microresonators (cf. Section 1.9). 
Based on our simulations we developed an open-source graphical user interface 
using MATLAB [83] that allows numerical simulation of Kerr-comb and soliton 
generation in microresonators. 


1.8 Characterization of DKS in crystalline mi- 
croresonators 


The ability to stably tune into soliton states as described in the Section 1.6 


allows for characterization of the solitons’ spectral and temporal properties in 
crystalline microresonators MgF 2 . Precisely as predicted by theory, the single 
soliton state has a sec/i-squared spectral envelope from which a pulse dura¬ 
tion in the femto-second regime can be inferred. The latter can be confirmed 
via direct temporal characterization. Here, frequency resolved optical gating 
using a nonlinear Michelson-interferometer fSHG-FROGl [HI] is used, cf. Fig¬ 
ure s This method also allows for the reconstruction of the complex pulse 
enveloDe|85j. The high pulse repetition rate inherent to microresonators allows 
the FROG interferometer delay to extend over several pulse round-trip times, 
i.e. consecutive pulses can be cross-correlated. 
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Figure 1.11: Numerical simulations of dissipative Kerr soliton formation in a 
crystalline microresonator, (a) Intracavity power (corresponding to the trans¬ 
mission signal in Figure when mirrored horizontally) during a simulated 
laser scan (101 simulated modes) over a resonance in a microresonator. The step 
features are clearly visible. The light gray lines trace out all possible states of 
the system during the scan. The unshaded area corresponds to the area where 
temporal DKS can exist, the light shaded area allows for breather solitons with 
a time variable, oscillating envelope; no solitons can exist in the dark shaded 
area. b. Optical spectrum and intracavity intensity for different detuning values 
(1-5) in the laser scan. 
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Figure 1.12: Characteristics of a temporal dissipative Kerr soliton in a crys¬ 
talline resonator. (a) Optical spectrum showing the characteristic sech^- 
envelope. (b) Magnified part of the spectrum, resolving the individual comb 
lines of which it is composed, (c) The low noise radio-frequency beatnote at 
14.09 GHz corresponds to the comb line spacing and the soliton pulse repetition 
rate, (d) SHG-FROG trace revealing femto-second pulse duration. The pulse 
to pulse separation of 71 ps corresponds to the pulse repetition rate of 14.09 
GHz and the pulse duration of 194 fs can be inferred (in agreement with the 
spectral width). (RBW: resolution bandwidth, GF: center frequency) 


Besides stable single soliton states, also stable multi-soliton states can be 
generate using the laser tuning method. Figure 1.13 compares the optical spec¬ 
tra and SHG-FROG traces of a single and two different multi-soliton states. A 
salient feature of multi-soliton states is the modulated spectral envelope arising 
from the interference of the individual soliton. Based on the spectral envelope 
the relative positions of the solitons can be reconstructed (cf. Section 1.11). 
Once a multi-soliton state is generated inside the cavity, the envelope pattern 
does not change over several hours, implying that the temporal separation be¬ 
tween the solitons in the microresonator does not change on this time scale. 



Figure 1.13: Spectral and temporal characterization of multi-soliton states in 
a crystalline resonator, (a) Optical spectra of a single soliton (top) and two 
multi-soliton states (middle, bottom), (b) SHG-FROG traces corresponding to 
(a). 
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1.9 


Resonator mode structure and soliton 
formation 


The generation of bright temporal DKS requires a focusing nonlinearity and 
an anomalous group velocity dispersion (GVD), i.e. a linear growth of the res¬ 
onator’s FSR with the mode number. In this case the growth of the FSR can 
be compensated by nonlinear mode pulling, i.e. nonlinear optical shifts of the 
resonance frequency. In principle most microresonators fulfill this requirement 
in the (near-) infrared spectral regime. As opposed to optical fiber systems, mi- 
croresonators, however, can exhibit a rich and complex optical mode structure. 
Some modes are affected by higher order dispersion or are even characterized by 
normal GVD. Moreover, fabrication defects and asymmetries can lead to modal 
coupling and distortions of the mode spectrum due to avoided mode crossings. 
While temporal dissipative Kerr solitons are a single-mode phenomenon, the 
structure of the respective mode-family can be altered by the presence of the 
other modes. The understanding of the impact and influence of the mode- 
structure is a decisive step in understanding DKS in optical microresonators. 
In Ref. |35j this is studied experimentally and with the help of the numerical 
tools described in Section [T?7| First, the mode structure of the microresonator 
is precisely measured using frequency comb assisted diode laser spectroscopy 
|86) . Figure 1.14 shows a fraction of the recorded transmission spectrum of a 
crystalline resonator, revealing several mode families. The full data set is rep¬ 
resented in Figure [1.15^ . Individual mode families and their dispersion can be 
extracted from the full data set. An example of a mode family with almost per¬ 
fect anomalous GVD is shown in Figure [T.15| 3, top. In contrast. Figure [T.15| 3, 
bottom shows a another mode family where the globally anomalous GVD is 
locally modified by two avoided mode crossing. While in the first case solitons 
can be generated, this is not possible in the latter. The hypothesis that avoided 
mode crossing (if they occur in the spectral vicinity of the pumped mode) inhibit 
soliton formation is confirmed by numerical simulation |35j . Moreover, strong 
higher order dispersion can prevent soliton formation. Importantly, moderate 
contribution of higher order dispersion and avoided mode crossings further away 
from the pump wavelength do not inhibit soliton formation but manifest them¬ 
selves in the optical spectrum as qualitatively illustrated in Figure [ lT6 | Indeed, 
similar structures can be observed in the experimentally generated soliton spec¬ 
tra e.g. in Figure [T.12[ 
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Figure 1.14: Transmission spectrum of a MgF 2 crystalline microresonator with a 
FSR of approximately 14.09 GHz. The upward transmission spikes (values > 1) 
result from cavity-ringdown. Frequency comb assisted diode laser spectroscopy 
ensures precise calibration of the laser detuning (MHz level). 
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Figure 1.15: Mode structure of a MgF 2 resonator with a FSR of 14.09 GHz. 
(a) Two-dimensional Echelle-type representation where for all measured mode 
families the deviation u)^ — loq — ^ ■ Ili,avg of the resonance frequency from 
an equidistant Hi^avg -spaced frequency grid (Dl,avg is an approximate average 
FSR of all modes) is shown (plus some offset) in function of the mode number ji. 
Dots forming a continuous line represent a particular mode family. Different free 
spectral ranges correspond to different slopes of the lines, whereas dispersion and 
variation of the FSR show as curvature and bending of the lines. The dispersion 
can be strongly affected by mode crossings, (b) Two mode families have been 
extracted from the data set shown in (a). The upper one is characterized by an 
anomalous dispersion, the lower one exhibits two avoided mode crossings that 
induce deviations from the anomalous dispersion. 
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Figure 1.16: Mode-structure and spectral envelope, (a) Typical sech-squared 
envelope for the case of weak (a) and strong anomalous dispersion (b). (c) 

Higher order dispersion (such as non zero II3 ) leads to an asymmetric spectrum, 
dispersive wave emission (cf. Section 1.11) and a shift of the spectral soliton 
peak intensity away from the pump laser (soliton recoil), (d) Avoided mode 
crossings manifest themselves in a characteristic spectrally local variation of the 
spectrum. 


1.10 Using dissipative Kerr solitons to count the 
cycles of light 

The observation of Kerr solitons in microresonators, and the associated genera¬ 
tion of ultra-short femto-second optical pulses, provides an opportunity to count 
the cycles of light using an optical microresonator, i.e. to link phase coherently 
an optical laser frequency to the radio frequency domain. In order to achieve 
such a radio to optical frequency link[87] self-referencing has to be employed; a 
process where the comb’s offset frequency /o (cf. Figure O) is measured. To 
this end, an optical spectrum of sufficient spectral span is required that spans 
an either a full octave when using the 1/-2/ self-referencing technique or 2/3 
of an octave for the 2/-3/ self-referencing method [3H1E71IT3] (cf. Figure [lT 7| . 
While the soliton spectrum generated in a crystalline resonator is not sufficiently 
broad by itself, the short pulse duration enables (after power amplification in an 
erbium-doped fiber amplifier) spectral broadening [88] in highly-nonlinear opti¬ 
cal fiber. Figure [T. 18^ shows a soliton spectrum that has been broadened in this 
way [HHl- Here, 2/-3/ self-referencing is implemented allowing for detection of 


the soliton comb’s offset frequency /o (cf. 1.18i). Together, with the comb line 
spacing (i.e. the pulse repetition rate) /rep that can readily measured via direct 
photo-detection (cf. |1.18| : ), the frequencies constituting the soliton comb are 
fully defined. In this way, the temporal dissipative Kerr soliton allows linking 
optical frequencies (e.g. the pump laser wavelength or any other comb line) to 
the radio (or microwave) frequencies /o and /rep. 
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Figure 1.17: Illustration of 2/-3/ self-referencing. If the frequency comb spec¬ 
trum spans more than two thirds of an octave, the second and third harmonics 
of blue and red wings of the spectrum can be overlapped. The difference fre¬ 
quency beatnote between the two harmonics yields the comb’s offset frequency 
/o. 
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Figure 1.18: Counting the cycles of light by self-referencing a soliton based 
frequency comb, (a) Soliton spectrum and nonlinearly broadened spectrum. 
The broadened spectrum spans more than two thirds of an octave and allows 
for 2/ — 3/ self-referencing, (b) Magnified part of the broadened spectrum. 
The line spacing is the same as for the soliton spectrum (14.09 GHz), (c) 
Pulse repetition rate beatnote measured via direct photo-detection of the comb 
spectrum, (d) Offset frequency signal measured via a modified 2/ — 3/ self- 
referencing technique, where two transfer laser are used for signal enhancement 
(see Ref. [89] for details). (RBW: resolution bandwidth, CF: center frequency) 
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1.11 Dissipative Kerr solitons and soliton induced 
Cherenkov radiation in a Si 3 N 4 photonic chip 


The discovery of temporal DKS in crystalline resonators, as described in the 
previous chapters not only enables generating coherent optical frequency combs 
in a reliable manner, but also opened the ability to use time domain broaden¬ 
ing methods to achieve a coherent link from the RF (or microwave) domain to 
the optical frequency domain. A key emerging research question has been if 
solitons can also be generated in other microresonator platforms. One particu¬ 
larly promising platform for applications are Si 3 N 4 integrated microresonators. 
This platform allows for the integration of coupling waveguide and resonators 
on a microchip. Indeed, it was demonstrated that by carefully optimizing the 
fabrication process, the influence of detrimental avoided mode crossings can be 
strongly reduced, enabling the generation of DKS in Si 3 N 4 microresonators [36] . 
Here, as in the case of crystalline resonators, the soliton formation is indicated 
by discrete steps in the resonator transmission. Similar to the method described 
above in detail in Section[T 6 |a suitable laser scan allows for stable and thermally 


locked solitons in Si 3 N 4 microresonators. Figure 1.19 shows the generation of a 
single soliton state in a Si 3 N 4 microresonator. The bandwidth of the generated 
spectrum is much wider than what has been observed in crystalline resonators. 
As an immidiate consequence higher order dispersion terms become relevant; in¬ 
deed, the third order dispersion causes the Kerr soliton to emit what is known as 
Soliton-Cherenkov radiation |51| (or dispersive wave), i.e. a radiative tail in the 
time domain. It is important to note that the process can also be understood in 
the frequency domain as successive FWM [90]. Dispersive wave formation has 
been predicted particularly for Si 3 N 4 resonators, where for a coherent spectrum 
a spectrally sharp dispersive wave is expectedjnT], as it is indeed the case (cf. 
Figure 1.19). For a resonator with third order dispersion (i.e. non-zero D 3 ), 
this occurs approximately near the mode number /r « — 3^ (three times the 
distance of the zero dispersion point from the pump wavelength, /x « — ^) [35] . 


It is important to note that unlike the case of a soliton propagating along 
an optical fibers where the dispersive wave is generally incoherent with the 
soliton (i.e. propagating at a velocity different from the soliton’s velocity), the 
coherence can be maintained in case of the microresonator owing to the periodic 
boundary conditions. This is similar to the coherence of the dispersive wave in 
supercontinua that are generated by a periodic train of ultra-short pulses. At 
the same time, the presence of a dispersive wave can also occur in the case of 
a high noise, non soliton state. The presence of a dispersive wave therefore is 
not a proof of a coherent Kerr soliton state [31]. Note, that the coherence of 
the spectrum in Figure 1.19 has been verified by recording the soliton pulse 
repetition rate beatnote at 190 GHz as well as the heterodyne beatnotes of the 
soliton spectrum with external CW lasers at various wavelengths (including the 
dispersive wave spectral region). Based on the spectral width a soliton pulse 
duration of approximately 30 fs can be estimated. 





29 


wavelength (nm) 



-200 -100 0 100 200 
rel. mode number jj 


Figure 1.19: Single Kerr soliton generation in a Si 3 N 4 microresonator. The 
upper panel shows the single soliton spectrum that covers a spectral bandwidth 
of 2/3 of an octave. The black line in the background is the spectral envelope 
obtained by the simulation described in Section |1.7[ The lower panel shows 
the measured dispersion (dots) free of avoided mode crossing and the dispersion 
over the full spectral span as obtained through hnite element simulation. The 
dispersive wave forms at the wavelength of approximately 1.9 micron where the 
phase matching condition — wq — /i ■ Di « 0 is fulhlled. The zero dispersion 
point (ZDP) is marked by a vertical, dashed line. The inset in the upper panel 
show a heterodyne beatnote between the dispersive wave and an external laser. 
Its narrow width proves the coherence of the dispersive wave. 
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Figure 1.20: Generation of multi-soliton states and soliton induced Cherenkov 
radiation in a Si 3 N 4 microresonator. The characteristic spectral modulations 
occur due to the interference of solitons at different positions in the microres¬ 
onators as indicated in the insets. 

It is also possible to generate multiple solitons in the Si 3 N 4 microresonators. 
As in the case of the crystalline resonators, once a multi-soliton state is gener¬ 
ated, the modulation of the spectral envelope, resulting from the interference of 
the solitons, remains stable for hours. The spectral envelope function /(/i) of 
this interference is given by 


/(m) = 


i=i 


where (pj corresponds to the relative angular position of the j**' soliton. From the 
envelope of the spectrum the single soliton spectrum and the relative positions 
of the solitons can be reconstructed as shown in Figure [T.20| 


1.12 Summary 

The observation and stable generation of temporal DKS in optical microres¬ 
onators provides a major impetus to the field of microresonator frequency combs. 
It provides a long sought method for the generation of widely spaced frequency 
combs, with smooth spectral envelope and with low phase noise. In addition. 
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soliton induced Cherenkov radiation provides a method to increase the band¬ 
width and transfer the coherence of the comb into the normal GVD regime. 
Together with the ability to accurately predict the comb performance via nu¬ 
merical simulations, this implies that the synthesis of octave spanning frequency 
comb spectra directly from a CW laser using a microresonator is a possibility. 
Solitons in microresonators may make frequency comb metrology ubiquitous by 
enabling compact comb generators that exhibit mode spacings in the 10 — 100 
GHz regime. Potential applications include frequency metrology (i.e. frequency 
measurements), data transfer, microwave signal generation, spectroscopy, opti¬ 
cal sampling and arbitrary optical waveform generation. 
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